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Based on the recent data in NUBASE2012, an improved empirical formula for evaluating the 
a-decay half-lives is presented, in which the hindrance effect resulted from the change of the ground 
state spins and parities of parent and daughter nuclei is included, together with a new correction 
factor for nuclei near the shell closures. The calculated a-decay half-lives are found to be in better 
agreements with the experimental data, and the corresponding root-mean-square (rms) deviation is 
reduced to 0.433 when the experimental Q-values are employed. Furthermore, the Q-values derived 
from different nuclear mass models are used to predict a-decay half-lives with this improved formula. 

It is found that the calculated half-lives are very sensitive to the Q-values. Remarkably, when mass 
predictions are improved with the radial basis function (RBF), the resulting rms deviations can 
be significantly reduced. With the mass prediction from the latest version of Weizsacker-Skyrme 
(WS4) model, the rms deviation of a-decay half-lives with respect to the known data falls to 0.697. 

PACS numbers: 23.60.-|-e, 21.10.Dr 


I. INTRODUCTION 

Alpha-decay is a very important process in the field 
of nuclear physics. It was first observed as an unknown 
radiation by Becquerel in 1896 and further formulated 
empirically by Geiger and Nuttall in 1911 [l|. After¬ 
wards, applying the theory of quantum mechanics in the 
field of nuclear physics, Gamow 0 and Condon and Gur¬ 
ney i independently described the spontaneous a decay 
as a quantum tunnelling effect through the potential bar¬ 
rier leading from the parent nucleus to the two emitted 
fragments: the a particle and the daughter nucleus. As a 
powerful tool, a decay can be used to investigate the low- 
energy structure of unstable nuclei, such as the ground- 
state energw the ground-state half-life, the shell effects, 
and so on pl-ll^. On the experimental side, the obser¬ 
vation of a-decay chains from unknown parent nuclei to 
known nuclei has been a reliable method used to identify 
different superheavy elements (SHEs) and isomeric states 
as well |13l - llq . 

Up to now, on the basis of Gamow’s theory, the ab¬ 
solute a-decay width has been estimated by many theo¬ 
retical calculations [13, [H, [2ll - l^ , which employ various 
approaches, such as the cluster model [iTl - f^ . the dei rsity 
dependent M3Y (DDM3Y) effective interaction [13, , 

the generalized liquid drop model (GLDM) [13]. Differ¬ 
ent from the cluster model, also some other theoretical 
models have been proposed in the pursuit of a micro¬ 
scopic description of a decay, such as the shell model 
and the fission-like model (26l - l^ . 

Moreoyer, many simple empirical formulas were also 
exploited to analyze the a decay [lol - l^ . Among these 
studies, Viola and Seaborg proposed a semi-empirical for¬ 
mula to analyze a-decay half-liyes for the heayy elements 
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{A > 140) 1^. which has been often used up to the 
present day |37l. Id^j. In 2000, Royer deyeloped other 
simple analytical formula by fitting on a completed set 
of 373 a emitters [s^j, which can well reproduce the ex¬ 
perimental half-lives of the favored a-decay. Based on 
the above two works, Sobiczewski and Parkhomenko [s^j 
also presented a simple phenomenological formula for de¬ 
scribing a-decay half-lives of heavy (above ^°®Pb) and 
superheavy nuclei. In previous works, the ground state 
spins and parities of parent and daughter nuclei were con¬ 
ventionally ignored. However, when the spin and parity 
values of parent and daughter nuclei are different, then 
the emitted a particle carries out nonzero angular mo¬ 
mentum 1. Because of this, the orbital moment of emit¬ 
ted a particle should be taken into account in an ac¬ 
curate approach for a decay. In this way, Denisov and 
Khudenko [s^j explored a carefully updated and selected 
partial a-decay half-life data set of 344 ground-state-to- 
ground-state (g.s.-to-g.s.) a transitions, and presented 
sets of simple relations for evaluating the half-lives of a 
transitions. After that, Royer [s^ proposed other sets 
of analytical formulas for logj^QTi /2 depending or not on 
the angular momentum of the a particle from an adjust¬ 
ment on the similar experimental data set mentioned in 
Ref. [sl] . Recently, Dong et al. [llj extended the Royer’s 
formula by taking into account the contribution of the 
centrifugal barrier, and then proposed a novel law, named 
the improved Royer’s formula, for calculatin g a -decay 
half-lives, which was proved to work well [^,[43]. Very 
recently, Ren et al. [ijl proposed a new Geiger-Nuttall 
law where the effects of the quantum numbers of a-core 
relative motion as well as the possible effect of angular 
momentum and parity of a particle are naturally embed¬ 
ded in the law. 

In the empirical formulas, the alpha-decay half-life is 
determined by the proton number Z, the number of nu¬ 
cleons in nucleus A, and the value of decay energy Q, 
which can be derived from the nuclear mass Q. Half- 
life calculations are very sensitive to the choice of Q- 
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values, so reliable theoretical predictions of the nuclear 
mass, leading to precise Q-values, are essential to study 
a transitions when these nuclear masses can not be de¬ 
termined experimentally. In the past decades, many 
theoretical calculations were performed to extrapolate 
nuclear masses. One conventional method is the local 
mass relations, which have a high precision of predic¬ 
tion for nearby nuclei, such as the Garvey-Kelson rela¬ 
tions [i^ , residual proton-neutron interactions [13, > 

Coulomb-energy displacement [13, [131 , and systematics 
of a-decay energies [5l|. The other method relies on 
global mass models, which are usually believed to have 
a better ability of mass extrapolation for nuclei far from 
the known region, such as the macroscopic-microscopic 
finite-range droplet model (FRDM) [13] , the Weizsacker- 
Skyrme (WS) model [131 , the microscopic Hartree-Fock- 
Bogoliubov (HFB) theory with a Skyrme force [13] and 
the relativistic mean-field (RMF) model [s^ j. 

In this paper, we synthetically consider the centrifugal 
effect and the hindrance of a emission with odd values 
of I, and then renovate the improved Royer’s formula 
by introducing two new empirical terms. The effect of 
angular momentum and parity of a particle is embed¬ 
ded remarkably, together with a new correction factor 
for nuclei near the shell closures. Comparisons about 
the root-mean-square (rms) deviation are performed, and 
the numerical results show that our formula can well re¬ 
produce the experimental half-lives of 341 a transitions 
between the ground states of nuclei. Moreover, we also 
consider different mass models to provide the Q-values 
for a-decay half-lives. Especially, when the radial ba¬ 
sis function (RBF) approach is canonically introduced to 
improve the precision of the mass models. A decrease of 
the rms deviations is obtained, which indicates that the 
calculated a-decay half-lives after embedding the RBF 
approach agree better with the experimental data. 

This paper is organized as follows. In section II, the 
theoretical method of calculating a-decay half-lives is 
briefly described. In section III, the numerical results 
of a-decay half-lives as well as some detailed discussions 
are given, including some comparisons between different 
empirical formulas and different mass models with (with¬ 
out) the RBF approach. Finally, we render a concise 
summary in section IV. 


state of parent nuclei and various states of daughter nu¬ 
clei for correct extracting of the half-lives for the g.s. to 
g.s. a decay. 

According to the updated information of Audi 
et al. [53], there are totally about 700 nuclei for all possi¬ 
ble a transitions with different i?-values. The nuclei with 
well-defined experimental i?-values are adopted in this 
paper. As a result, our data set contains 341 nuclei with 
their i?-values being in the region of 1% < R < 100%. 
These nuclei can be divided into four categories: 123 even 
Z-eveii N (e-e), 93 even Z-odd N (e-o), 79 odd Z-even 
N (o-e), and 46 odd Z-odd N (o-o) nuclei. 


B. Spin and parity selection rnle 


In general, for e-e case, the spin and parity values of 
parent and daughter nuclei are usually ignored. However, 
for the other three cases, the transitions may occur with 
different spins and parities of the parent and daughter 
nuclei and, consequently, the a particle may take away 
an nonzero angular momentum 1. Accordingly, the effect 
of the orbital moment of emitted a particle should be 
taken into account in an accurate approach for a decay. 
We select I in accordance with the spin-parity selection 
rule, similar as that in Refs. [H, HH], 


I 


min 


A, 

A, 

A, 

A, 


for even Aj 
1 for even Aj 
for odd Aj 
1 for odd Aj 


and TTp = TTd, 
and TTp y^TTd, 
and TTp^TTd, 
and TTp = TTd, 


( 1 ) 


where Aj = \jp—jd\, jp, T^p, jd and tt^ are spin and parity 
values of the parent and daughter nuclei, respectively. It 
is worth noting that the orbital angular momentum of 
the emitted a particle can have several values according 
to the selection rule. However, for the sake of simplicity, 
in the following calculations, the angular momentum I 
of the emitted a-particle is endowed with a minimum 
value ^min- In addition, the values of spin and parity of 
parent and d aug hter nuclei in this paper are taken from 
Refs. Bllllli- 


C. Analytical formula for the a decay half-lives 


II. THE THEORETICAL METHOD 
A. Input experimental data 

In this paper, a recent a-decay data of 341 nuclei is 
taken from NUBASE2012 [^, and then a fitting proce¬ 
dure on this data set leads to an improved empirical for¬ 
mula. Moreover, the following study is initially restricted 
to the g.s. to g.s. a transitions, which are known as par¬ 
tial a decay of the parent nucleus. In view of this, it is 
necessary to take into account experimental values of the 
branching ratio (i?) for a transitions between the ground 


In 2000, by fitting the experimental lifetimes of 373 
emitters with i?-values close to 100%, Royer gave an el¬ 
ementary analytical formula [33| . 

r7 

^ogif^{Tii2) = a + bA^^^yfz + ■ ( 2 ) 

Here a, b and c are the fitting parameters, A and Z are 
the mass number and charge of parent nucleus, respec¬ 
tively, Ti /2 is the half-life of a decay and given in sec¬ 
onds while Q the corresponding decay energy and given 
in MeV. This formula can well reproduce the experimen¬ 
tal half-lives of the favored a-decay. However, it dose not 
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TABLE I: The parameters of our empirical formula found for 
e-e, e-o, o-e, and o-o ranges of nuclei, respectively. 




a 

b 

c 

d 

e- 

■e 

-25.432 

-1.146 

1.577 

— 

e- 

■o 

-26.591 

-1.171 

1.639 

1.123 

o-e 

-27.747 

-1.093 

1.620 

0.829 

o- 

-0 

-28.460 

-0.984 

1.573 

0.970 


consider the angular momentum in both unfavored and 
hindered a. transitions. Just because of this, the above 
formula was reconsidered by Dong et al., and converted 
into a new version in terms of the centrifugal contribution 
to the a-nucleus potential [s^, 


logio(’^i/2) — o. + + ^jQ- 

\/(A-4)(Z-2)A-2/3' 


(3) 


In this paper, after detailed investigation about the 
explicit dependence of the a-decay half-lives, we propose 
a different formula, 


logio(Ti/2) = a + bA^/^y/Z + 


cZ 


+ 


y/Qa 

V^(A-4)(Z-2)A-2/3 


+ 5. (4) 


The four parameters a, b, c and d can be obtained by fit¬ 
ting the new experimental data set, and then expressed 
in Tab. HI It is interesting that we find the coefficient 
is a value less than 1 for odd I in both o-e and 
0-0 cases, so this may imply that the fourth term in the 
formula proposed in Eq. (3) may overestimate the cen¬ 
trifugal effects for these two cases. 

In Eq. (4), the first three terms are the same as those 
in Refs. [321138|| . the fourth term is similar to the corre¬ 
sponding one in Ref. , which can be attributed to the 
distinct contribution of the centrifugal potential 
to the total a-nucleus potential at small distances be¬ 
tween daughter nucleus and a particle while I ^ 0 [s^. 
The difference is that we have introduced an additional 
factor into this term to account for the hindrance 

of the transition with the change of parity. It is such 
an exponential form that the formula can well reflect 
the changes of alpha-decay half-lives originated from 1. 
Due to the rearrangement of single-particle orbits and of 
the nuclear spin orientation, the structure hindrance is 
relatively weak and, consequently, will not be dis¬ 
cussed here. The last empirical term is a phenomeno¬ 
logical correction factor and defined as: S = 0.5 for 
49 < Z < 51, 81 < Z < 83, 49 < Af < 51, 81 < iV < 83 
or 125 < N < 127. S is introduced to mock up the sys¬ 
tematic deviations of calculated a-decay half-lives with 
respect to the experimental data near shell closures. 


In addition, the rms deviation of the decimal logarithm 
of the a-decay half-life, in this paper, is defined as 




N 


N 


5^[logio(T,7--) - logio(7’r/T )]"- (5) 


N is the number of nuclei used for evaluation of the rms 
deviation. The Q-value for the a-decay half-life can be 
calculated using experimental mass data as Q = Mp — 
{M(i Ma), Mp, Md, and are the masses of parent 
and daughter nuclei and a particle, respectively. 


III. RESULT AND DISCUSSION 



FIG. 1: (Color online) Logarithms of the ratios between 
theoretical a-decay half-lives calculated with different formu¬ 
las and experimental ones versus the mass number A of the 
parent nucleus. The red diamonds denote the results of this 
work, the blue circles and the olive squares correspond to the 
results calculated with two formulas proposed by Dong et al. 
and Royer, respectively. 


As proposed above, Royer’s formula works well in the 
favored a-decay [13 • However, when it is used to eval¬ 
uate the half-lives for nuclei with Z and N crossing the 
Z = 82 OT N = 126 shell closures, respectively, a dra¬ 
matic large deviation can be observed between calculated 
half-lives and experimental ones, which is drawn clearly 
in Eig.[T]with the squares corresponding to the results of 
Royer’s formula in Eq.(2). It can be seen in Eig. |T]that 
all squares are far away from the dotted line, the cal¬ 
culated values of a-decay half-lives are always less than 
the experimental data for both Z = 81 — 83 isotopes 
and N = 125 — 127 isotones. There are about a half of 
the whole isotopes and isotones with the ratios between 
calculated values and experimental ones being beyond a 
factor of ten (log^g 10 = 1). The systematic deviation is 
very apparent. 

Then after including the centrifugal effect, Royer’s for¬ 
mula was extended to the unfavored a-decay with a 
new form proposed in Eq.(3) [^ . With this formula. 
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the same calculations are performed on both the iso¬ 
topic chain oi Z = 81 — 83 and the isotonic chain of 
N = 125 — 127, and we plotted the corresponding results 
with the filled circles in Fig.[TJ It is obvious that the val¬ 
ues of logiQ(T^^|° /T®Ji^* ) for the most nuclei (denoted 
by filled circles) are in the range from —1 to 0 and the 
whole points land close to the dotted line. It means the 
calculated a-decay half-lives with Eq. (3) agree better 
than those deduced from Eq. (2) with the experimental 
data. However, the deviation is still systematic with the 
mean absolute value of log;^g(r*y 2 °'/’^i/^*) being about 
0.5. Accordingly, the systematic behavior of deviation is 
decreased to some extend but not completely redressed. 
Then how to overcome the deviation error between the 
calculated half-lives and the experimental data? 

As we know, accurate calculations of a transitions 
should take into account the spins and parities of par¬ 
ent and daughter nuclei and the angular momentum of 
the emitted a particle. Particularly for the transition 
with large change of spin and with the change of parity, 
the factor can well simulate the hindrance of a 

emission. Moreover, when Z {N) goes across the shell 
closure at Z = 82 (A^ = 126), the effect of the closed 
shell results in a decrease of the a-preformation factor, 
which has been shown in the experimental analysis [^ . 
that may be the reason of the existence of the large de¬ 
viations proposed above. In this way, we introduce two 
new empirical terms to remedy the systematic deviation, 
and a new analytical formula is proposed in Eq. (4). 
Applying our formula to both Z = 81 — 83 isotopes and 
N = 125 — 127 isotones leads to a desired effect, which 
can be seen in Eig. [T] with the red points. One can see 
that most of the red points float around the dotted line, 
the absolute values of logig(T)^y 2 °'/^i/ 2 *) generally 
less than 0.4 for almost all of the parent nucleus. Accord¬ 
ingly, the systematic behavior of deviation is decreased 
to a large extent. 

For further insight of the dependence on I and S', a 
comparison between the experimental a-decay half-lives 
and the calculated ones for whole 341 nuclei with the 
three formulas proposed above is performed, and is drawn 
in Fig. in which the results of the empirical formulas 
are calculated with the parameter values from the present 
work. In Fig. [51 Panel (a) denotes the decimal logarithm 
of Ty ^°-for 341 nuclei with Royer’s formula pro¬ 
posed in Eq. (2). One can see that around the shell 
closures {Z = 82 and N = 126), the calculated half- 
lives are much shorter than the experimental ones, the 
ratios between calculated values and experimental ones 
are beyond a factor of ten (logig(T^^|° /T®J’ 2 *’) = ~1)- 
The rms deviation for the full data set is 0.659. In panel 
(b), the formula proposed in Eq. (3) is applied. It is 
obvious that the agreements with the experimental data 
are visible around the shell closures. The amounts of 
discrepancy also decrease in other fields, and the rms 
deviation is reduced to 0.487, all of which indicate the 
extending of Royer’s formula by considering the centrifu- 
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Neutron number 

FIG. 2: (Color online) Decimal logarithms of T‘^^' 
for 341 nuclei derived from the different formulas, for which 
calculations are performed with the present parameter set of 
this work. Panel (a) denotes the case with Royer’s formula 
proposed in Eq. (2), panel (b) corresponds to the case with 
the improved Royer’s formula proposed in Eq. (3), and panel 
(c) is designated the result obtained with our formula pro¬ 
posed in Eq. (4), respectively. Dotted lines denote the magic 
numbers. 


gal effect is successful. Panel (c) is designated the result 
obtained with our formula proposed in Eq. (4). In this 
panel, it is worthwhile to note that systematic behavior 
of deviation is rapidly decreased. The absolute values 
of log;^Q(T*^ 2 °'/^i/:r*’) is optimized in many fields, and 
averagely about 0.27 for nuclei near the shell closures. 
Our rms deviation is only 0.433. Hence, the introduction 
of the term S and parity correction improves 

strongly the efficiency of the formula proposed in Eq. 
(4), even though they are semi-empirical. By the way, 
at the top right-hand corner in Panel (a)-(c), the ratios 
between calculated values and experimental ones are be- 
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TABLE II: Comparisons between the rms deviations of the 
experimental and the calculated a-decay half-lives for the full 
data set as well as for e-e, e-o, o-e, and o-o subsets, respec¬ 
tively. The first column contains the symbols for correspond¬ 
ing formulas. 


Total e-e e-o o-e o-o 


{N = 341) {N = 123) {N = 93) {N = 79) {N = 46) 


Eq. (2) 

0.587 

0.298 

0.713 

0.610 

0.810 

Eq. (3) 

0.481 

0.298 

0.574 

0.486 

0.637 

Eq. (4) 

0.433 

0.267 

0.521 

0.434 

0.574 

Ref. [3^ 

0.536 

0.350 

0.489 

0.539 

0.689 

Ref. [37] 

0.561 

0.285 

0.509 

0.548 

0.747 


yond a factor of ten for some isotopes or isotones and not 
decreased after introducing I and S. The reason may be 
that the SHE experiments are very difficult and usually 
few decay events are observed. Hence the experimental 
error bar is relatively large in the measurement of both 
decay energies and half-lives. 

The detailed numerical results for the full data set are 
listed in Tab. m as well as for e-e, e-o, o-e, and o-o sub¬ 
sets, respectively. Note that in order to study the effect 
of S and I on the accuracy of the formula, the parame¬ 
ters used in either Royer’s formula (HI or the improved 
Royer’s formula [s^ have been recalculated by fitting 
the new experimental data set. After this, it is found 
that the corresponding rms deviations deduced from the 
above two formulas are 0.587 and 0.481, less than the 
old a of 0.646 and 0.504 with the parameters taken from 
Refs. [13, [13, respectively. It means that the accuracies 
of these two formulas are improved apparently, but still 
lower than that of our formula with the standard cr-values 
being only 0.433. Moreover, a further comparison is per¬ 
formed between our formula with the well-known empir¬ 
ical formulas presented by Denisov and Khudenko |35| 
and Royer in 2010 [131, respectively. The calculations 
are also performed with the new parameters determined 
by fitting the updated experimental data set, and then 
lead to the corresponding cr-values of 0.536 and 0.561, 
respectively. The detailed rms deviations for both full 
data set and four subsets are displayed in Tab. |TT]in the 
fourth and fifth lines, from which a similar conclusion is 
attained in comparison with our formula. That is our 
analytical formula has the smallest values of the rms de¬ 
viations for both full data set and most of the subsets 
except for the e-o case, in which the rms deviations re¬ 
sulted from Refs, [la, 1131 are smaller. All these indicate a 
higher accuracy in our formula than either one in Ref. [^ 
and Ref. [33 |. 

To examine the predictive ability of nuclear mass 
model on a-decay half-life, the improved formula ex¬ 
pressed in Eq. (4) is employed to calculate a-decay half- 
lives with the theoretical Q-values. In this work, eight nu¬ 
clear mass models are taken into account, i.e., the FRDM 
model [13 , the latest version of WS (WS4) model [111 , 



Neutron number 


FIG. 3: (Color online) The same as Fig. (2] but for different 
four mass models without (with) the RBF approach. The 
left four panels (al-dl) correspond to the case without the 
RBF approach, while the right ones to the case with the RBF 
approach. 


TABLE III: The same as Tab. [Ill but for different calculated 
mass tables. 

Model Model-I-RBF 


FRDM 

1.622 

0.968 

WS4 

1.163 

0.697 

HFB27 

1.617 

1.128 

RMF 

5.822 

1.511 

Bhagwat 

1.184 

0.862 

DZIO 

1.644 

0.979 

DZ31 

1.597 

1.035 

KUTY 

1.184 

0.890 


the recent version of the HFB model (HFB-27) 113 . the 
RMF model with TMA effective interaction [53, the 
Koura-Tachibana-Uno-Yamada (KTUY) model [61|, the 
Duflo-Zuker formulas [H, [ill and the Bhagwat formula 
[Hi- For simplicity, the deviations of log 2 o(T‘^ 2 °'/’^i/ 2 *) 
for four mass models are drawn, and without loss of gen¬ 
erality, FRDM, WS4, HFB-27 and RMF are taken as ex¬ 
amples (see Fig. [3]). The corresponding numerical results 
of rms errors for eight mass models are listed in Tab. IHII 
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As can be seen in Tab. imi the descriptions of a-decay 
half-lives for these mass models are generally within the 
same order of magnitude except for the RMF model. The 
smallest rms deviation of 1.163 is obtained from the WS4 
model, which indicates the WS4 model has the best ac¬ 
curacy of Qa prediction among mass models considered 
here. However, the RMF model significantly overesti¬ 
mates the a-decay half-lives for nuclei around N = 126 
and in some other fields (see panel (dl) in Fig. [3]), so it is 
necessary to improve the mass precision of RMF model 
for reliably predicting a-decay half-lives. 

An efficient method to enhance the predictive power 
of mass models is the radial basis function (RBF) ap- 
proach, which was hrst introduced into this field by Wang 
et al. [6^. The basic formulas of the RBF approach have 
been detailed in our previous works 1^ . So we will 
not repeat it anymore hereafter. For convenience, the 
mass model improved by the RBF approach is denoted 
with the Model-I-RBF henceforth, e.g. FRDM-I-RBF and 
RMF-I-RBF. After employing the RBF approach, the 
new a-decay half-lives can be deduced with the optimized 
Q values derived by the above mass equation, and sub¬ 
sequent comparisons between the new a-decay half-lives 
and the experimental data are performed and exhibited 
in Tab, mil One can see from Tab. lIIIl that the systematic 
rms deviations decrease remarkably after using the RBF 
approach for all the eight mass models. Taking RMF 
model for an instance, the rms deviation is 5.822 before 
applying the RBF approach, and then significantly de¬ 
creases to 1.511 with refinement by the RBF approach. 
From Fig. [21 one can also see that the calculated a-decay 
half-lives with RBF approach result in better agreements 
with the experimental ones for RMF model than that 
without RBF approach. The effect of systematic cor¬ 
rection is impressively clear. Meanwhile, there are the 
common cases for the other mass models, all of which 
provide strong support for the reliability of the RBF ap¬ 
proach and its usefulness for eliminating the discrepancy. 


IV. SUMMARY 

In summary, by htting the new experimental data of 
341 nuclei, we have proposed an empirical formula for 
systematically calculating the half-lives of a transitions 
between the ground states of parent and daughter nuclei. 
Within this new formula, we consider the hindrance 
effect resulted from the changes of spin and parity of 
parent and daughter nuclei and a phenomenological 
correction factor for the nuclei near the shell closure. 
The calculations of a-decay half-lives with some other 
well-known empirical formulas are also performed. By 
comparing, it is found that our formula results in the 
smallest values of the rms deviation for both full data 
set and any subset. In other words, the precision in our 
formula is higher than that in the previous methods. 
Furthermore, we apply the improved formula to eight 
calculated mass tables. Using the Q-values resulted from 
these mass tables, the subsequent half-life outcomes 
indicate that WS4 model has the best accuracy of Qa 
prediction. Then in order to improve the predictive 
ability of the mass model, the RBF approach is adopted. 
As a result, a rapid decrease of the rms deviation 
is obtained, which indicates the calculated a-decay 
half-lives, after embedding the RBF approach, are in 
better agreements with the experimental data. 
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